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ABSTRACT 

There  is  the  largest  concentration  of  forest  resources  in  India  along  with  equally  large  concentration  of  tribal  and 
rural  population.  Tribal  population  and  rural  population  are  highly  dependent  on  forest  resources  for  their  livelihood. 
Degradation  and  depletion  of  forest  resources,  increase  poverty  and  suffering  among  the  rural  population,  therefore  it  is 
imperative  to  educate  the  tribal  population  to  use  the  forest  products  efficiently  and  make  efforts  for  sustainable 
management  of  forest  resource  biomass.  Proper  awareness  of  tribal  population  is  demand  of  the  time.  Keeping  these 
objectives  in  mind,  we  have  formulated  a mathematical  model  consisting  of  biomass,  tribal  population,  trained  tribal 
population,  awareness  programme  to  aware  and  train  the  tribal  population  for  efficient  use  of  forest  products.  Local  and 
global  stability  analysis  of  the  mathematical  model  along  with  the  persistence  of  the  system  is  checked  using  the  theory  of 
nonlinear  ordinary  differential  equations.  Analytical  results  obtained  are  justified  numerically  through  numerical 
simulation.  The  important  parameters  are  investigated  and  a variety  of  variables  to  change  in  these  parameters  is 
determined. 

KEYWORDS:  Forestry  Biomass,  Population,  Stability,  Persistence 

1.  INTRODUCTION 

We  are  living  in  the  twenty  first  century.  The  world  has  seen  spectacular  political,  social,  cultural,  economic  and 
scientific  progress  in  this  century.  But  when  we  look  at  cost  paid  by  us  to  achieve  this  progress  it  seems  that  what  we  have 
gained  is  much  less  than  the  loss  faced  by  the  environment  due  to  its  over  exploitation.  Over  the  years,  there  has  been 
progressive  pressure  on  the  environment  and  the  natural  resources,  the  alarming  consequences  of  which  are  becoming 
evident  in  the  form  of  frequent  natural  disasters  in  the  form  of  earthquakes,  droughts  and  floods  throughout  the  world. 
These  consequences  have  direct  impact  on  the  lives  of  the  poor  who  are  directly  dependent  on  natural  resources. 
Demographic  and  economic  growth  has  historically  driven  deforestation,  forest  exploitation,  and  agricultural  demand 
(South,  1999;  Wright  and  Muller-Landau,  2006;  Cohen,  2014a,  b),  and  such  growth  will  continue  in  the  developing  world. 
Our  forest  wealth  is  deteriorating  due  to  over-grazing,  over-exploitation  both  for  commercial  and  household  needs, 
encroachments  by  the  tribal  and  rural  population.  Economic  growth  in  the  developing  world  is  projected  to  double  global 
consumption  of  forest  products  by  2030  (WWF  and  IIASA,  2012:  Ch4.  p.  9).  Other  unsustainable  practices  include  shifting 
cultivation  and  other  developmental  activities  such  as  roads,  buildings,  irrigation  and  power  projects  in  the  rural  areas  as 
well. 
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On  one  side  there  are  deforestation  and  an  increase  in  unsustainable  practices  in  rural  areas,  the  other  side  includes 
an  ever  increasing  movement  of  poor  families  to  towns.  In  fact,  Lack  of  opportunities  for  employment  in  villages  has  led  to 
an  ever  increasing  movement  of  resource-poor  families  in  the  town  leading  in  expansion  of  urban  slum  areas.  Illiteracy  and 
child  labour  are  persisting  in  urban  areas.  Records  say  that  the  global  rate  of  forest  loss  has  decreased  since  2010  to  3.3 
million  hectares  (Mha)  or  0.08%  annually  (Keenan  et  al.,  2015).  Similarly,  rates  of  afforestation  are  steady  or  rising  not 
just  in  temperate  countries,  were  planted  forests  have  long  been  integral  elements  of  the  forest  estate,  but  also  in  the 
tropics,  where  the  extent  of  planted  forests  has  nearly  doubled  since  1990  (Payn  et  al.,  2015).  Notwithstanding  this 
progress,  forest  conversion  for  agriculture  has  not  reduced  and  may  lead  to  eventual  forest  loss  (Van  Lierop  and  Lindquist, 
2015). 

Hence,  seeing  this  scenario,  awareness  or  educating  the  masses,  especially  in  the  tribal  and  rural  population  is  one 
of  the  main  tools  that  will  help  to  perform  transition  to  a greener  Earth  and  saving  the  environment  from  continuous 
deterioration.  Awareness  helps  in  disseminating  innovations  in  technology  to  the  people  who  need  it  in  actual.  It  is  the 
ability  of  the  society  to  generate  new  ideas  and  approaches  to  solve  environmental  related  problems.  The  demand  for  new 
skills  is  growing  and  in  order  to  accomplish  this  task,  the  systems  of  education  and  professional  training  will  need  to 
provide  well-trained  and  highly  qualified  personnel. 

Keeping  all  these  things  in  mind,  we  formulate  a nonlinear  mathematical  model  using  ordinary  differential 
equations,  considering  different  equations  for  resource  biomass,  tribal  population,  aware  or  educated  tribal  population  and 
awareness  programme.  There  have  been  many  studies  investigating  the  impact  of  industrialization  and  increasing 
population  pressure  on  the  environment.  Many  researchers  have  investigated  the  depletion  of  resource  biomass  by 
overgrowing  population,  toxicants  and  industrialization.  Shukla  et  al  (1989,  1996,  and  1997)  has  studied  the  effects  of 
population  and  pollution  on  the  depletion  of  forestry  resources.  Dubey  and  Dass  (1999)  have  proposed  and  analyzed  a 
mathematical  model  to  study  the  survival  of  species  dependent  on  a resource,  which  is  depleted  due  to  industrialization. 
Dubey  and  Narayanan  (2010)  studied  a mathematical  model  to  demonstrate  the  effect  of  industrialization,  population,  and 
pollution  on  the  depletion  of  a renewable  resource.  M.  Agarwal  et  al.  (2010)  proposed  a ratio  dependent  mathematical 
model  on  the  depletion  of  forestry  resource  biomass  due  to  industrialization  pressure.  Shukla  et  al.  (2003)  studied  the 
effects  of  primary  and  secondary  toxicants  on  the  resource  biomass.  Shukla  et  al.  (2011)  studied  the  effect  of  technology  on 
the  conservation  of  forestry  resource  biomass.  It  was  further  observed  that  in  the  tropics,  where  most  forest  change  has 
occurred,  deforestation  due  to  smallholder  agricultural  expansion  has  given  way  to  large-scale,  enterprise-driven  forest 
conversion  (Hecht,  2005;  Rudel,  2007;  Asner  et  al.,  2013).  Sayer  and  Cassman  (2013)  argued  that  future  agricultural 
demand  for  land  could  be  met  without  significant  forest  loss  if  agricultural  innovations  were  deployed.  Sloan  and  Sayer 
(2015)  gave  a review  of  Forest  Resources  Assessment  of  2015  showing  positive  global  trends,  but  still  it  was  observed  that 
forest  loss  and  degradation  persist  in  poor  tropical  countries.  Going  through  the  above  literature,  it  was  observed  that  none 
of  them  considered  tribal  and  trained  tribal  population  along  with  the  effect  of  awareness  on  them  for  sustainable 
management  of  our  resources,  which  is  novelty  our  work.  We  aim  at  finding  the  impact  of  training  tribal  population  of 
conservation  of  resource  biomass.  Our  paper  is  classified  as  follows:  In  section  one  brief  introduction  about  the  problem  is 
given,  in  section  two  mathematical  models  are  given.  Section  three  includes  boundedness  of  the  system,  section  four;  five 
and  six  illustrate  equilibrium,  stability  (both  local  and  global)  analysis  and  persistence  of  the  model.  In  section  seven  and 
eight  numerical  simulations  and  conclusion  of  the  model  respectively  is  given,  followed  by  the  list  of  references. 
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2.  MATHEMATICAL  MODEL 


Let  B is  biomass  density  of  forest  resources,  P is  tribal  population,  T is  aware  tribal  population  and  A is  an 
awareness  programme  to  train  the  tribal  population.  It  is  assumed  that  resource  biomass  are  growing  in  intrinsic  growth 

rate  r,  tribal  population  is  consuming  the  resourc.iomass  at  the  rate  & . A)  is  the  rate  at  which  resource  biomass  is 
benefited  by  the  awareness  program  launched  among  the  tribal  population.  The  intrinsic  growth  rate  of  tribal  population  is 

s-  Carrying  capacity  of  resource  biomass  and  tribal  population  is  given  by  K and  respectively.  The  rate  at  which 

tribal  population  gets  awareness  about  sustainable  consumption  of  forest  resources  is  Y ,s  0 and  Y()  are  the  migration  rate 
of  tribal  and  aware  tribal  population  to  the  urban  areas  in  search  of  better  employment  opportunities.  V is  the  rate  at  which 
awareness  programme  is  launched,  Awareness  programme  is  launched  proportional  to  the  tribal  or  rural  population  of  the 

system.  V()  is  the  failure  rate  of  awareness  program. 
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Where  B( 0)  > 0,  P( 0)  > 0, 7(0)  > 0,  A(o)  > 0. 


3.  BOUNDEDNESS  OF  THE  SYSTEM 


(2.1) 


In  analogy  to  the  population  dynamics,  it  is  very  important  to  observe  the  consequences  that  restrict  the  growth  of 
the  population.  In  this  sense,  the  study  of  boundedness  of  the  solution  of  system  around  different  steady  states  is  very  much 
needed.  For  this,  we  find  boundedness  on  the  system  in  the  following  lemma: 

Lemma  3.1:  The  set  O = {(£,  P.T,A):0<  B < K,0<  P < Pm  ,0  < T < Tm  and  0 < A < Am  }, 


Where 


p _{s-ul+aK)Kl  y PmAm 

s ’ " jU2-maK 


and  A 

m 


_ nP 
7 ’ 


with  conditions  s + aK  > > tnaK 

is  the  region  of  attraction  for  all  solutions  initiating  in  the  interior  of  the  positive  octant. 

Proof:  Let  ( B , P,  T,  A ) be  solution  with  positive  initial  values  (Bq  , P{] , T{) , A0  ) .Define  a function 
From  the  first  equation  of  the  system  (2.1) 
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dt  I K 


According  to  comparison  principle,  it  follows  that 


B = K. 


From  the  second  equation  of  the  system  (2.1),  we  get 


dP  f P 

— <aP1 +aKP-yPA-/JlP, 

dt  K, 


dP  / „ sP2 

— <(s+a  K-fijP--— 
at  K. 


According  to  comparison  principle,  we  have 


{s+a  K-jUt  )K1 


with  condition 


s + a K>  Hi- 


From  the  forth  equation  of  the  system  (2.1),  we  find  it 


-pA 


According  to  comparison  principle,  we  get 


From  the  third  equation  of  the  system  (2.1),  we  get 


— < maTK  + y P A -u,T, 

i d / m m / 2 7 

dt 


-p^  rP,  Am-(jU2-maK)T, 
dt 

T ?PmK 

m Hi  ~ maK  ’ 


with  condition 
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H2  > maK. 

Therefore,  all  solutions  of  the  system  (2.1)  enter  into  the  region 

a = {{B,P,T,A) : 0 < B < K,0<  P < Pm,0  < T < Tm  andO  < A<  Am}, 


where 


p _{s-Jul+aK)Kl  T = 


yPmAm 
jU2  - maK 


and  A = 


rjPm 

7, 


With  conditions  s + CXK  > //j , /A  > maK 

is  the  region  of  attraction  for  all  solutions  initiating  in  the  interior  of  the  positive  octant. 

This  completes  the  proof  of  the  lemma. 

In  the  next  Section  we  present  the  equilibrium  analysis  of  the  model. 

4.  EQUILIBRIUM  ANALYSIS 

The  system  (2.1)  has  four  nonnegative  equilibria  namely;  (o,  P,  T,  A)and  E3  {B*,  P*,  T*,  A *).The 
existing  of  equilibrium  points  E0  (0,0,0,0)  and  El  ( A",0,0,0),  is  obvious.  We  show  the  existence  of  other  equilibria  as 
follows 

The  existence  of  En  (o  ,p,t,a) 


Here  P,T  and  A are  the  positive  solutions  of  the  system  of  algebraic  equations  given  below 


f 

s 
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i-A 

K. 
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-y  A-  Ei 


0, 


(4.1) 


yPA-jUj  = 0, 

rjP-vA-Tj  0f  = 0. 

From  the  equations  (4.1)  and  (4.2),  we  get 

-={s-£±-yA)KL 

s 

with  condition 

S > //,  + y A. 

yKl{s-nl-yA)A 

Mi 


(4.2) 

(4.3) 

(4.4) 

(4.5) 

(4.6) 
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Using  the  values  of  P and  T from  the  equations  (4.4),  (4.6)  with  condition  (4.5)  in  equation  (4.1),  we  get 
j_A2±2]a22-4A1A3 


Where 

A =srj0r2Kl,A2  = yrjKljul+julsrjl  + r/0yKxs{s  -//,),  A3  =//,^l(.v -//,). 

This  completes  the  existence  of  E2 . 

The  existence  of  E3  ( B P*,  T*,  A *) 

Here  B*,P*,T*  and  A * are  the  positive  solutions  of  the  system  of  algebraic  equations  given  below 


r 1 -aP*  -maT*  = 0, 

l K 


s 1 + aB*-y  A*-Jul  = 0, 

K, 


maT  * B * +y  P * A*  -ju2T*  = 0, 
tjP  * -rj]  A * -r/0T*  = 0. 


From  equation  (4.9),  we  get 


P*  =—  1-  — 

a K 


- maT * = f(B\Tt). 


From  equation  (4.10),  we  have 


1 i St"  i 1 / \ 

A*  = \saKl  - sr  + B*  — + Kxa 2 ~ii3Kxa  + sma2T*\  = fAB*,T* ). 

K.ayi  \K  ) J 


Using  the  equations  (4.13)  and  (4.14)  in  the  equations  (4.1 1)  and  (4.12),  we  have  following  equations 


v ' K \K 


+ Kial  j|  + T‘{-  JulK1a+  rsma- ma(sK,  a - sr  - jUt  Kt  a)} 


+ B*<—  — +Kia1  \ + —(sKla-sr-julK1a)[-—\—+Kla2  S*2 -m2aisT’3  + 
J aK  \ aK\K  J 

-(sKfl-sr-  //,  K{a) . 


H2(b\T*)=TH  1-— r/maT  — ^-{sKia-sr+B'[—+Kla2\-nlKla 

# ^ K J K^cxy  \K  J 


2r*}-7  of. 


From  equation  (4.15)  we  note  the  following 
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When  T — > 0,  then  B — > Ba  where 

cX2  + c2b:  + c3  = 0. 


(4.17) 


Where 


r ( sr  ^ 


C,  = — — + K,a2  ,C,=KC , (sK,a - sr-ju,K,a), C3  = K2C,  - KC2 

aKyK  ) aK 


Thus  equation  (4.17)  has  unique  positive  solution  if  following  conditions  hold 
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Now  from  equations  (4.19)  and  (4.20),  we  get 

dH, 

dB * _ dT” 


dT * dH,  ' 

dB * 


(4.21) 


dB  n 

It  is  clear  that  > U,  if  either 

dT 


dH,  dH, 

— — r > 0 and  f < 0,  or 

dB  dT 


(4.22) 


dH,  dH, 

—4  < 0 and  —4  > 0 

dB  dT 


hold 


From  equation  (3.16)  we  get  following 
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When  7 — > 0,  then  B — > Bb  where 

, _ //  rKKj+  ?]l  Ksr  — r/i  KKta(\  + //, ) 
rj\sr+  KK^^+rj  rKj 

With  condition 

? ) rKKx  y+7)lKsr>7)l  KK1  a(  1 + juJ. 


377  2 


—ijma  - 


Tjxsma 2 
Kxay 


-%■ 


377,  rjr  r/l  ( sr  2 

— f = — + K,a 

dB  aK  Kxay\K 


Now  from  equations  (4.25)  and  (4.26),  we  get 

dH2 

dB * = dT* 
dT * 377,  ' 

3 


It  is  clear  that 


377, 

3/L 


> 0 and 


dB 
dT * 

377, 


3r* 


> 0,  if  either 


< 0,  ( 


(4.23) 


(4.24) 

(4.25) 


(4.26) 


(4.27) 


(4.28) 


^<0and^f>0 


3 B 


3 T 


hold. 

( * m*  i 

B ,T  ) if  in  addition 

to  the  conditions  (4.22)  and  (4.28),  the  inequality  7>fl  < Bh  holds.  Knowing  the  value  of  B and  7 , the  value  of  r and 
A can  be  calculated  from  equations  (4.13)  and  (4.14).  This  completes  the  existence  of  E3  (B*,  P*,T*,  A *)  . 


5.  LOCAL  STABILITY 

To  discuss  the  local  stability  of  the  system  (2.1),  we  compute  the  variational  matrix  of  the  system  (2.1).  The 
entries  of  the  general  vocational  matrix  are  given  by  differentiating  the  right  side  of  the  system  (2.1)  with  respect  to 
B,  P,  T and  i.e. 
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v{e)= 


^11  ^12  ^13  0 


t2l  t22 


0 L 


*31  *32  *33  *34 

0 t42  t43  tu 


Where 

2rB  2 sP 

tn  = r aP  - mod , tl2  = -aB,  tl3  = -maB,  t2l  = aP,  t22  = s + aB  - yA  - //, , 

K K{ 

tM  = -yP,  t31  = mod , t31  = yA , t33  = maB  - //2 , t34  = yP,  t42  = 7,  t43  = -tj0  , t44  = -rjx . 

Accordingly,  the  linear  stability  analysis  about  the  equilibrium  points  , i = 0,1, 2, 3 gives  the  following  results: 

1 . The  equilibrium  point  E0  is  unstable  manifold  in  B — P plane  and  stable  manifold  in  T — A plane. 

2.  The  equilibrium  point  E]  is  unstable  manifold  in  P — direction  and  stable  manifold  in  5 — T — A plane. 

3.  The  equilibrium  point  E2  is  unstable  manifold  in  B — direction. 

The  stability  behaviour  of  equilibrium  point  E3  is  not  obvious.  However,  in  the  following  theorem  we  give 
sufficient  conditions  for  equilibrium  points  E3  to  be  locally  asymptotically  stable.  For  proof,  see  Appendix  A. 

Theorem  (5.1):  Let  the  following  inequalities  hold: 

(y  A*  f < 2S  — [Un  ~ maB " ),  (5.1) 

K, 


(P'-riXfc^K. 


(5.2) 


Then  an  equilibrium  point  E3  is  locally  asymptotically  stable. 

Global  Stability 

The  following  theorem  characterizes  the  global  stability  behavior  of  equilibrium  point  E3  .For  proof,  see 
Appendix  B. 

Theorem  (5.2)  Let  the  following  inequalities  hold: 

rn2a 2 (rj0Tm  - P’vf  < 2r71^F  fa  - maB * ),  (5.3) 

K 

f K < ^r—  {ju2  - maB * ).  (5.4) 

KJJ  o 

Then  an  equilibrium  point  E3  is  globally  stable  in  the  region  D. . 


www.iaset.us 


editor@iaset.us 


122 


Archana  Singh  Bhadauria,  Rachana  Pathak  & Manju  Agarwal 


6.  PERSISTENCE 

Theorem  (6.1):  Assume  that  r > aPm  + maTm , s + aBmm  > yAm  +Hi>Hi  > m a 5min  and  VP^n  > Vo  Tm  ■ 
Here  K,P,T  and  A are  upper  bounds  of  the  populations  B,  P,  T,  A respectively  and  always  positive.  Then  the 
system  (1)  becomes  persists. 

Proof:  From  the  first  equation  of  the  system  ( 1 ),  we  have 


dB  / x rB 

— >{r-aPm-  m a Tm  )B — . 

dt  K 

According  to  lemma  (3.1)  and  comparison  principle,  it  follows  that 

= K(r  - a Pni  - rnaTm) 

^min 


(6.1) 


(6.2) 


With  condition  r > OC  P + mC(Tm  , 5min  remains  always  positive. 
From  the  second  equation  of  the  system  (1)  and  equation  (6.2),  we  have 

^ > (5  + aBmm-yAm  - //,  )P - P2 . 
dt  Kt 

According  to  lemma  (3.1)  and  comparison  principle,  it  follows  that 
p _ {s  + aBmn-yAm-^)Kx 

min 


(6.3) 


(6.4) 


With  condition s + C(Bmin>  yAm  + //, , Pmm  remains  always  positive. 
From  the  last  equation  of  the  system  (1 ) equation  (6.4),  we  have 

dA 

— ^ VPmm  -VoTm-Vi  A. 

dt 

Using  the  comparison  principle,  it  follows  that 

iVP^n-V  0Tm) 


A - = 


Vi 


(6.5) 


(6.6) 


With  condition  T)  P mm  > Vo^m  Amn  remains  always  positive. 

From  the  third  equation  of  the  system  (1)  and  equations  (6.2),  (6.4),  (6.6),  we  get 

^7  ^ rP„Am„  ~ (M2  ~ (6-7) 

dt 

According  to  lemma  (3.1)  and  comparison  principle,  it  follows  that 
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7’ ■ = 


- m a Bn 


(6.8) 


With  condition  /LL-,  > niCXBmm  , 7’m|n  remains  always  positive. 

This  completes  the  proof  of  the  theorem.  Then  the  system  (1.1)  becomes  persists. 

7.  NUMERICAL  SIMULATION 

The  global  stability  of  the  non-linear  model  system  (2.1),  in  the  positive  octant,  is  investigated  numerically. 
The  numerical  integration  of  model  systems  (2.1)  is  carried  out  for  various  choices  of  biologically  feasible  parameter 
values  and  for  different  sets  of  initial  conditions.  In  all  the  cases  being  considered  here  the  data  are  chosen  such  that  the 
persistence  conditions  are  satisfied.  Model  system  (2.1)  is  solved  numerically  using  the  Runge-Kutta  method.  The 
dynamic  behavior  and  its  corresponding  time  series  of  the  model  system  (2.1)  are  decided  on  the  following  data  set: 

r = 10,  K = 150,  a = 0.02,  m = 0.001, .?  = 3,  Kt  = 100,  y = 0.001,  ^ = 0.001,  n2  = 0.3, 77  = 0.03, 

7,=  0.04, 7/0  =0.02.  (7.1) 

The  interior  equilibrium  point  of  the  model  system  (2.1)  corresponding  to  the  feasible  parameter  values  as  (7.1)  is: 

£3(100.80,  163.95,  53.05,  96.43). 

Figures  1,  2,  3 and  4 are  drawn  to  study  important  parameters  in  the  system.  Figure  1 is  plotted  to  display  the 
variation  of  resource  biomass,  tribal  population,  aware  tribal  population  and  awareness  programme  with  time.  It  manifests 
the  local  stability  of  the  system  at  an  interior  equilibrium  point  as  each  population  approach  the  equilibrium  value  as  time 
tends  to  infinity. 


Figure  1 
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Figure  4 
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Figure  5 

DISCUSSIONS  AND  CONCLUSIONS 

Tribal  population  and  rural  population  are  highly  dependent  on  forest  resources  for  their  livelihood.  Degradation 
and  depletion  of  forest  resources,  increase  poverty  and  suffering  among  rural  population,  therefore  it  is  imperative  to 
educate  the  tribal  population  to  use  the  forest  products  efficiently  and  make  efforts  for  sustainable  management  of  forest 
resource  biomass.  Proper  awareness  of  tribal  population  is  demand  of  the  time.  Keeping  these  things  in  mind,  we  propose 
and  solve  a mathematical  model  to  study  the  impact  of  awareness  program  among  tribal  and  rural  population.  Our  study 
includes  model  formulation  followed  by  qualitative  analysis  of  the  formulated  nonlinear  mathematical  model.  Qualitative 
analysis  of  the  model  includes  Equilibrium  analysis,  local  and  global  stability  analysis  and  persistence  of  the  system. 
Numerical  simulation  is  performed  to  justify  the  analytical  findings  and  graphs  are  plotted  to  study  the  variation  of 
important  variables  of  the  system  with  time  for  different  parameters.  Analytically  conditions  for  the  local  and  global 
stability  of  the  system  are  determined  and  these  conditions  are  further  justified  numeric  by  considering  a hypothetical  set  of 
parameter  values.  Local  and  global  stability  of  the  system  is  displayed  through  the  graph  also.  Through  the  graph  it  is 
observed  that  as  the  rate  of  launch  of  awareness  program  in  the  system  increases  the  equilibrium  value  of  resource  biomass 
increases.  This  implies  that  as  people  get  awareness  about  sustainable  management  of  resources,  exploitation  of  resources 
decreases  and  hence  the  original  equilibrium  level  of  resource  biomass  increases.  Thus,  awareness  program  helps  to 
increase  the  equilibrium  level  of  resource  biomass  and  it  is  imperative  to  be  motivated  both  by  the  government  and  public. 
In  addition,  it  is  observed  that  the  density  of  resource  biomass  decreases  gradually  by  continuous  use  of  resource  biomass 
without  any  awareness  program.  However,  when  awareness  program  is  launched  then  resource  biomass  retain  its 
equilibrium  value  and  approaches  the  same  equilibrium  level  as  if  it  is  untouched  by  the  tribal  population. 
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Appendix  A.  Proof  of  Theorem  (5.1) 

To  prove  Theorem,  we  first  linearize  model  (1)  by  substituting 

B-B  +b.P=P  + p,T  = T*  +t,A  = A*  + a. 


■ ■ “ a=0,m=0,Y=0 
a=0.02,m=0,Y=0 
~ a=0.02, m=0.001, y=l 
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Where  b,  p,  T,  a are  small  perturbations  around  an  equilibrium  point  E , . We  get  a following  linearize  system 

db  rB* 


dt 


K 


-b-  aB  p-maB  z. 


dp  sP * 

— = p + aP  b—yP  a, 

dt  Kt 


(A.l) 


— = maT*b+  y A* p+  y P*a  + (maB*  - jU0)t, 
dt 


da 

dt 


= ?lp-Tila-Tj<>T. 


Then  we  consider  the  following  positive  definite  function: 


W = —b 1 + — p2  +—T1  +— a". 
2 2 2 2 


Differentiating  W with  respect  to  time  t along  the  linearize  system  (A.l),  we  get 

(PAJ  . 2 i 7 2 2 2 

= —wnb  + wl2bp  + wnbr+  w^4ap  + w22pr  +w34ar  - w22p  - w22t  - wua~ . 

dt 


(A.2) 


(A.3) 


Where 


rB 

W11  =—  ^12 


Sit  P* 

■ a(n j P - B ),  w13  = ma(n2T  - B ),  w22  = — - — , w33  = n2  (jd2  - maB  ), 

K , 


iv. 


44 


= nfh,  w24  = n37/-yntP\  w34  = n2yP  -n3?j0,  w23  = yn2A*. 


ist  u • B />’  yPB 

Now  choosing  n = — n = — n = ' 

P ' T tl0T 


(A.4) 


Using  (A.4),  equation  (A.3)  reduce  in  following  form 


dAV  . 2 2 2 2 

= -wnb  + w24ap  + w2ipf-  w12p  -w33f  - w44a  . 

dt 


(A.5) 


dW 

Sufficient  conditions  for to  be  negative  definite  are  that  the  following  inequalities  hold: 

dt 


(r  A*  )2  < [p2  ~ maB* ), 

Ki 


(A.6) 


KtfoyT 


(A.7) 


This  completes  proof  of  Theorem  (5.1 ). 
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Appendix  B.  Proof  of  Theorem  (5.2) 

For  finding  the  condition  of  global  stability  as  E 3 we  construct  the  Lyapunov  function 


H = B-B* -B*ln^-  + c 
B 


r-r  -r  in 


+ - 


(B.l) 


Differentiating  H with  respect  to  time  t along  the  solutions  of  model  (1),  we  get 

jg  = (g-g)jg  +c>-f)jP  (r_r)g 

dt  B dt  P dt  dt  dt 

Using  a system  of  equation  (1),  we  get  after  some  algebraic  manipulations  as 

d^  = -UB-B*)2 -SMP-P*)2 -c^-maB'lT-r)1  -c^A-A')2 

at  K 

+ {-  ma  + mac  2 TfB - B * \T-r)+(y  c2  P*  - c3rj0 \t -T*\a~  A*)+ 

(< c3rj  - c j y){A  - A*  \P  - P* )+  (Cla  - a )(b  - B*  \p  - P* ) ■ + y c2  a(t  - T*  \p  - P* ) . 

For  simplicity,  we  choose  c _ i c _ c _ Z 

1 ’ 2 3 77' 

Using  (B.4),  equation  (B.3)  can  be  written  in  following  form 

d^  = -^(B-B*f  -^(P-P*f  -c2(jU2-maB*lT-T*)2 -ctf^A-A*)2 
dt  K Aj 

+ {-  ma  + mac2T}{B - B*)(T  - T")+  yc2  A (t-T'\p~P*). 


(B.2) 


(B.3) 

(B.4) 


(B.5) 


dH 

A sufficient  condition  for to  be  negative  definite  are  that  the  following  inequalities  hold, 

dt 


m2a2(tjJm-  P tjf  < lrTJ^P  (j/i  - maP  ), 

K 


fK„<2yl^(ji2-maB"). 


(B.6) 

(B.7) 


Thus  the  model  (2.1)  is  globally  stable  in  the  region  D.  . 
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